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Abstract 

A refined hydrochemical model for single-bubble sonoluminescence (SBSL) is presented. The 
processes of water vapor evaporation and condensation, mass diffusion, and chemical reactions 
are taken into account. Numerical simulations of Xe-, Ar- and He- filled bubbles are carried out. 
The results show that the trapped water vapor in conjunction with its endothermic chemical 
reactions significantly reduces the temperature within the bubble so that the degrees of ionization 
are generally very low. The chemical radicals generated from water vapor are shown to play an 
increasingly important role in the light emission from Xe to He bubbles. Light spectra and pulses 
are then computed from an optically thin model. It is found that the resulting spectrum intensities 
are too small and the pulse widths are too short to fit to recent experimental results within stable 
SBSL range. Addition of a finite-size blackbody core to the optically thin model improves the 
fitting. Suggestions on how to reconcile the conflict are given. 
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I. INTRODUCTION 



The discovery that acoustic energy can be converted to hght through an oscillating air 
bubble trapped in water [l| has trwered extensive studies on single-bubble sonolumines- 
cence (SBSL) (see reviews in Refs. |2|,|3[). Under certain conditions, a stable and regular flash 
of blue- white light with a width of about 40—350 ps is emitted during the violent collapse 
of the bubble in each acoustic cycle 0,0,0]. Among many candidate mechanisms of SBSL 
light emission, the model that identified SBSL light emission as thermal bremsstrahlung 
and recombination radiation from an optically thin bubble seemed to predict the widths, 
shapes, and spectra of the emitted light fairly well under certain simplified hydrodynamic 



frameworks 



10|. Similar mode 



variable blackbody core 
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s included those which used finite opacity to define a 
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15| , whose suitability was refiected in the good fitting 
. However, different approximations in modelling 



by a smaller-than-bubble blackbody 
the physical-chemical processes, and particularly, uncertainties under the extreme conditions 
inside a sonoluminescing bubble, may give diverse predictions that plague the validity of a 
light emission model. Therefore, more realistic hydrodynamic-chemical modelling and more 
critical tests of the light emission models under the refined hydro-chemical framework are 
necessary. 

Past studies delineated the important effects of diffusive transport, surface tension, anc. 



compressibility of the surrounding liquid on SBSL 



Earlier studies 
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considered the influence of evaporation-condensation phenomena on the bubble dynamics. 



Sochard et al. j22| and Gong et al. 



23l | coupled the bubble dynamics with the water vapor 



dissociations. Recently, the study of sonochemsitry was extended from lower temperature 
to higher temperature situations where SBSL emerges 0, [2^. Yasui 



24l | presented a 



model of SBSL that accounts for evaporation-condensation process at the bubble interface 
and water vapor chemical reactions. It was later stressed by Storey and Szeri 2q] and 



others 
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that water vapor reduce the temperatures inside the SL bubble 
signiflcantly by reducing the compression heating of the mixture and through primarily 
endothermic chemical reactions. In some of these models, spatial uniformity of the bubble 
interior was assumed and inter-molecular mass diffusion was not properly accounted for. 
As a consequence, such models tended to underpredict the amount of trapped water vapor 
during the rapid collapse. 
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In the full hydrochemical numerical study by Storey and Szeri j26[, the consequences of 
water vapor inside strongly forced argon bubbles were investigated in detail. The interaction 
of nonlinearity of the volume oscillations, mass diffusion, and nonequilibrium phase change 
at the bubble wall resulted in excess water vapor trapped in the bubble during the violent 
collapse. The amount of trapped water vapor was more than that predicted by the simple 
model [2^. Akhatov et al. [3^ accounted for the occurrence of supercritical conditions 
of condensation and studied laser-induced cavitation bubbles. The effects of water vapor 
diffusion in different noble gas bubbles were studied by Xu et al. Isill . where shock waves were 
found to occur only in Xe bubbles. More recently, Toegel et al. j2S| studied the effects of the 
highly compressed conditions of SL bubbles on chemical equilibrium constants. They showed 
that high temperatures could be recovered due to the suppressed water vapor dissociations. 
In spite of the progress made, the direct consequence of sonochemistry on existing popular 
light emission models is to be revealed in a full hydrodynamic model. 

In this paper, we present a refined hydrodynamic model taking into account the chem- 
ical reactions and ionizations of the noble gas and water vapor mixture. The model is an 
extension of ou. p.eviou. one. QHQ. As done fo. a p.e a,-on bubble Q, tKe Nav.e. 



Stokes (NS) equations for the multispecies gas mixture in the bubble interior are coupled 
with a proper form of the Rayleigh-Plesset (RP) equation for the bubble wall, including 
the effects of liquid compressibility and heat transfer. The newly added feature is that the 
nonequilibrium processes of evaporation and condensation, species diffusion, noble gas ex- 
change between the bubble and the surrounding liquid, and dissociations of water vapor and 
ionizations of atomic species inside the bubble are all taken into account. The numerical 
scheme for solving the hydrodynamic equations is modified to be a semi-implicit one which 
allows for better numerical stability than original explicit scheme in Ref. [l8[. 

Detailed formulae in Eulerian framework are given. Numerical simulations are carried out 
for bubbles of He, Ar, or Xe gases. The effects of sonochemistry on current light-emitting 
models of SBSL N, MJ, are studied in detail through comparison with a calibrated 
experiment [3] • The main conclusions are: (i) the chemical reactions reduce the temperature 
within the bubble to such an extent that the degrees of ionization are generally very low; (ii) 
shock waves do not appear in He or Ar bubbles in the stable SBSL regime, but can occur 
in Xe bubbles only at higher driving pressures; (iii) chemical radicals generated from the 
water vapor contribute dominantly to the light emission of He bubbles; (iv) based on the 
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computed photon absorption coefficients, the hght spectra and pulse widths computed from 
the popular optically thin model can hardly be fitted to the experimental ones. Addition of 
a dynamic blackbody core to the optically thin model improves the fitting. Some suggestions 
on how to improve the optically thin model are outlined. 



II. HYDRODYNAMIC MODEL 



1511181 to mclude processes 



In this section, we extend our previous hydrodynamic model 
of evaporation and condensation on the wall and chemical reactions inside the bubble. The 
bubble is assumed to be spherically symmetric and is composed of the mixture of noble 
gas, water vapor, and reaction products. In addition to the NS and RP equations, the 
equations for the mass concentration of the dissolved noble gas and for the temperature of 
the surrounding water are also solved. 



A. Gas dynamics in the bubble 



1. The NS equations 



The bubble is assumed to contain A^-species gas mixture. Nonequilibrium chemical re- 
actions of the water vapor and ionizations of the monatomic species (Ar, H, and O) are 
considered. The maximum ionization level is taken as 3 for a noble gas, and 1 for H or O 
atom. For an Ar bubble, A^ = 13: Ar, Ar+, Ar2+, Ar3+, H2O, OH, H, H+, O, 0+, H2, O2, 
e~. The dynamics inside the bubble is described by the compressible NS equations, which 
represent the conservation of mass, momentum and energy. They can be written into a 
"conservative" form in the spherical coordinates 
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where pi is density of species i, p is the total density of gas mixture, P is the pressure, T is 
the temperature, u is the mass averaged velocity, e is the total internal energy of the mixture, 
Cj and hi are the internal energy and enthalpy of species i, respectively. Pi is the partial pres- 
sure, p and A are the viscosity and the thermal conductivity of the mixture, respectively, Yi 
and Xi are the mass fraction and mole fraction of species i, respectively. Mi is the molar mass 
of species i, M is the mean molar mass of the mixture, Dij is the binary diffusion coefficient 
between species i and j, Df^ is the mean diffusion coefficient of species i into the mixture, 
KJ is the thermal diffusion ratio, r^^ is the normal stress, Jj is the mass diffusion flux that 
must satisfy X^ili Ji = 0; cjj is the net mass production rate due to chemical reactions 
and ionizations that satisfies Y^.f^i oji = 0. Use of the mean diffusion coefficient is a practical 

nn 

approximation for computational efficiency [2a, However, to ensure global mass con- 
servation, a correction diffusion flux Jf = p, ^^I, + (X, -Yi)^ + Kf^] 



hi 



Ci + 
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dr I- dr 

is added to Jj in Eq.Q, as recommended by Ref. |32|- It can be easily shown that the 



modifled Jj satisfles Yld=i Ji — 0- 



2. Transport properties, equation of state and thermodynamic properties 

The individual transport properties (viscosity /ij, thermal conductivity Aj, thermal dif- 
fusion ratio kf , and binary diffusion coefficient Dij) are generally calculated based on 
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Chapman- Enskog theory [33|, |34, |35[. However, the transport properties of the gas mix- 
ture, n, A are determined by some empirical combination rules such as Wilke's semiempirical 
formula 34]. There is also difficulty in describing individual /i, and Dij for some reaction 
products due to lack of data. We shall let unavailable fii and Dij equal to other known 
ones, e.g., /ioh = /^H20;/^h = /^H2- Collision cross sections of ions, electron, and some radical 
species that are not available in Ref. 

dependent polynomial fitting js^. Once fii is known, Aj is obtained by modified Eucken 
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3J| can be determined using NASA temperature- 



model [33|, |3J|. The trouble lies in the determination of kf. In this regard, we only take 
into account the thermal diffusion between the noble gas and the water vapor, since ther- 
mal diffusion is only important in slow stage other than during collapse |26l. l35 1 . The high 
pressure corrections [34| are applied to /Xj and Aj similar to what Xu [3l| did. The transport 
properties of the mixture fi and A are obtained by using Wilke's semiempirical formula 
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The hydrodynamics of the bubble is affected by the equation of state. Here the gas 
mixture is modeled by a hardcore van der Waals equation of state that has the excluded 
volume but ignores the van der Waals force as the previous authors did 
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where R = J^iLi^i^i^^i = Ru/^u Ru is the universal gas constant, h = J^iLi'^i^i a 
simple combination of bi, with hi being the van der Waals excluded volume in m^/kg. The 
values of bi = biMi (in cm^/mol) are computed hj b = RTc/8Pc where and Pc are critical 
temperature and pressure of a species [31. When these critical parameters are not available, 
bi is taken as 4 times the spherical volume of the atomic or ionic radius. 

A well-posedness of equation closure requires P = /(e, pi, ■ ■ ■ ,p^), therefore, T must be 
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solved for from the given energy relation 
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where the internal energy is divided into the translational and rotational, the vibrational, and 
the ionization parts. The coefficient of specific heat at constant volume Cvi, is approximately 
assumed as follows |[24|: Cyi = | for monatomic gases such as Ar, Ar"*", H, O; Cvi = | for 
diatomic gases such as OH, H2, O2, and Cvi = | for other gases. e° is the reference energy, 
ef is the vibrational energy of molecules species, and ef is the ionization energy _of_ ion 
species. The reference energy e° is taken to be standard heat of formation at 298 K 



3. Chemical kinetics 

The chemical kinetics consists of the reaction mechanism and determines the net produc- 
tion rate of each species. For the chemical reactions of water vapor, we use the mechanism 
that were described in detail in Ref. js^. Only a subset consisting of eight elementary reac- 
tions is used, corresponding to the first eight ones used by Yasui '24^. The first 19 reactions 



of Ref. [2^ with additional species (HO2, H2O2) were also tried but the resulting temper- 
ature was found to be a little lower than that from the 8 reaction scheme. The processes 
of nonequilibrium collisional ionization and recombination are considered only for three 
monatomic species, noble gas, H, and O. The reason to choose H and O atoms is that they 
are quite ample in the water vapor dissociations, have lower ionization potentials, and can be 
treated using previous rate formulae. Ionizations of molecular species such as OH and H2O 
are believed to be more complicated, thus are not accounted for. The net mass production 
rate Uk due to chemical reactions of the water vapor is determined by the law of mass action 

= MkJ2Ki-Ki)Qi, (9) 

i=l 

AT / XI/' N y X u'.'. 

where is the net rate of progress of reaction i, and A^^ the total number of reactions. The 
forward and backward reaction rate constants for the ith reaction kf^ and kb- are given in 
Arrhenius form 

kj^ = A{T^( exp (-Cl /t) , = AlT^' exp (-C^/T) . (11) 



Since the above rate constants as given in Ref. |39| are based on ideal gas, the modification 
for a van der Waals gas as suggested by Toegel et al. is used. As the y h inted, we can 
derive the equihbrium constant based on fugacity for a van der Waals gas j34 1 
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Let the corrected forward rate frozen k'f, = kf^. We can obtain the 



corrected backward rate ki as 
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We will^ use k'j. and k'^^, in Eq. pOj) . The Effects of the above modification suggested in 



Ref. 



were found to suppress water vapor dissociation to some extent when compared 
with the raw rate constants kf^ and ki,.. 



When involving the third-body reaction Eq. (jlUj) becomes 
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where Z^^i is the third-body enhanced coefficient. Due to length limitation, the detailed 
formulas for net production rates and rate constants of collisional ionization, recombination, 
and three-body recombination are not given here. One can refer to Ref. [l^ |4y] for detail. 



4- Mass and heat exchange at the bubble wall 

The evaporation-condensation process and diffusion of the noble gas into the surrounding 
liquid are included. The net evaporation rate (mass per unit area and unit time) at the 
bubble wall is given as follows [2O, Q] '■ 

^ _ O^M ( -Psat(^f,int) _ ^Pv,mt j ^^^s^ 
v2^rR^ \ a/T; int \/TyAnt I 



where is the accommodation coefficient (evaluated using the formula in Ref. 2^) that 
shows which portion of water vapor molecules hitting the liquid surface is absorbed by this 
interface, is the gas constant of water vapor, Psat is the saturation vapor pressure at 
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liquid temperature T^^int, Pv,int is the partial pressure of water vapor, and T is the correction 
factor: 



r = exp(-fi2 
Q = 



'TX 



rfif. 



(16) 
(17) 



Pv,mt V 2 

In this study, the jump of temperature across the interface is assumed zero, thus T;^int = Tv;mt- 
Although Eq. (fT3j) is valid only below a critical point (for water ~ 647 K), it is used 
throughout the whole acoustic cycle for simplicity. The rate of mass diffusion of the noble 
gas dissolved in the liquid at the wall is 

dc 



rhd = AttR Di 



dr 
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where R is the bubble radius, Di is the diffusion coefficient of the noble gas, and c is the 
mass concentration of the noble gas dissolved in the liquid. 

The boundary conditions of species and energy at the bubble surface are derived by 
balancing the flux and source/ sink of an interface control volume with infinitesimal thickness. 
In doing so, we assume that there are no uptakes of radical species or surface reactions. The 
resulting boundary condition for gas species is 



Ji\R + mtotYi = rhefi + mafi 



where 

{1, i = water vapor 
0, i = other species 

The boundary condition of energy is 
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1, i = noble gas 
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where T/ and A; are the temperature and thermal conductivity of the liquid, respectively, L 
is the latent heat of evaporation of the liquid, and AH is the heat of solution of noble gas 



into the hquid |3J,|41|. 



The gas and liquid velocities at the bubble surface and the velocity of the bubble wall 
differ due to mass transfer. The boundary condition for gas and liquid velocities are 
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The mass and heat transfer boundary conditions are nonhnear. All the boundary conditions 
for the gas dynamics have to be coupled with the following motion, temperature, and noble 
gas concentration equations in the surrounding liquid. 



B. Motion, heat, and mass transport in the liquid 



The liquid flow outside the spherical bubble is accounted for with different approximations 
for motion and heat (or mass) transport, respectively. On the one hand, the Euler equations 
for the liquid motion can be reduced to the ordinary differential equation for the bubble 
radius known as the RP equation. The RP equation is coupled with the NS equations 
through stress equilibrium condition at the bubble wall. On the other hand, we assume 
that the fluid is incompressible when accounting for heat and mass transfer in the liquid. 
The separate treatment reduces the complexity of solving a fully coupled hydrodynamic 
equations at sacrifice that shock waves in the liquid can not be simulated well. 

Because the mass transfer at the bubble wall results in very small liquid velocity whose 
effect on the RP equation can be ignored, we let ui\r = R. Thus we can use a form of the 
RP equation that includes first order terms in the Mach number M = R/Cn, and 



allows for variable speed of sound in the water 
{1- M)RR+^(l- ^Mj R'^ = {1 + M) 



4|: 
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(23) 



Here subscripts b and 00 denote bubble wall and infinity, respectively, Ps(t) = —Pa sin(27r/t) 
is the pressure of the sound field with frequency / and amplitude Pa- For water, an equation 
of state of the modified Tait form 



P + B 



Pi 

Ploo 



(24) 



Poo + B 

is used with B = 3049.13 bar and n = 7.15. The enthalpy Hh and the speed of sound Cih of 
the liquid at the bubble surface are given by 
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b Plb 

The pressure Pib on the liquid side of the bubble surface is related to the pressure P{R,t) 
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on the gas side of the bubble surface by normal stress equilibrium condition 

P(R,t)-Trr\r=R = Plb + ^+'^ , (27) 

where t] is the dynamic viscosity and a is the surface tension. Their values depend on T^, as 



formulated in the Appendix of Refs. 



3G 



4l|, respectively. 

Both heat and mass transfer are taken into account although the former is found to be 
more important to the bubble dynamics. The equations for the water temperature and the 
mass concentration of dissolved noble gas take a similar form: 

dTi , dTi Xi d f_^dTi\ 



dt dr PiCpi r'^dr \ (9r / ' ^ 
dc dc ^ d f r,dc\ 

— +ui— = Di-—ir'—), (29) 



dt dr j.2Qj, \^ Qj, ^ 

where Cp^ is the specific heat at constant pressure of the liquid. In Eqs. and (j^ . the 
liquid velocity can be determined by the incompressible assumption 

= (30) 

The boundary condition for the water temperature are the continuity of heat flux Eq. ()2ip 
and Ti |^^^ = Too, and the boundary condition for the mass concentration is 

cU«=^^%^Pno(i?,t), cUoo=Coo, (31) 
-TO 

where cq is the saturated dissolved gas concentration at Too and Pq, and Pno is the partial 
pressure of the noble gas on the internal side of the bubble interface. 

C. Numerical method 

To exploit the advantage of a stationary Eulerian meshes, we use x = r/R{t) to transform 
the NS equations (0) to a form in fixed domain x G [0, 1] as done in our earlier work The 
transport equations for water temperature T/ and mass concentration c in domain r G [R, oo] 
are transformed into diffusion- type equations in domain z G [0, 1] through two consecutive 
coordinate transformations 45, 3| with the aid of Eq. (j3(Jj) . The details were given in 
Refs. Q,|42|. 

We apply the second-order upwind total-variational-diminishing (TVD) scheme j3l to 
the inviscid flux terms and the central difference to the diffusive terms of the NS equations. 
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The temporal discretization differs from Ref. |l8[ in tfiat we now use an Adams-Basliforth 
explicit scheme for the convective and spherical coordinate terms, and Crank-Nicolson im- 
plicit scheme for diffusive and chemical source terms. The implicit treatment is to overcome 
the stiffness problem due to diffusive transport and chemical source terms. The trapezoidal 
rule and central difference are used for the water temperature and concentration equations. 
A predictor-corrector method is used for the RP equation. We use 400 grid points for the 
NS equations and 100 points for the water temperature and gas concentration equations. 



III. OPTICAL POWER RADIATED BY THE BUBBLE 

n 

We compute SBSL based on the weakly ionized gas model of Hilgenfeldt et al. |8], which 
was most thoroughly studied and remarkably successful [lO.J. In more general cases with 
nonuniform bubble interior, the bubble has an optically thin radiating / absorbing outer shell, 
and may have a blackbody inner core when the opacity is large enough (lll.ll2|. We previously 
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applied this generic version to a pure argon bubble [IJ, In this paper, we apply the 
formulas of photon absorption coefficients 0, to the gas mixture in the bubble. The 
overall photon absorption coefficient is the sum of contributions from all species 

TV N 
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where is the absorption due to the free-free interaction of electron and ions, is 
the absorption due to free-free interactions of electrons and neutral atoms, and t^^^ is the 
absorption by bound-free ionization of already excited atoms. The ionization potentials 
used can be found in Ref. fs^. The bound-bound absorption is not accounted for and the 
modification to electron-neutral-atom bremsstrahlung Isil is not adopted for the time 
being. 

The finite opacity model given in Refs. jl^ . llS^ (without the G correction) is used to 
compute the total spectral radiance (power emitted per wavelength interval) of the bubble 
content at wavelength A 

Pl\t) = p levr/tf (r, t)Rl\r, t) exp {- 4°*(r', t)dr'^ r^dr (33) 



+ATTRlRl\R,,t)e^V [- j 
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where 

is the spectral emissive power (energy per unit time, wavelength interval, and projected 
surface area) with the Planck and Boltzmann constants h and ks, and the light speed in 
vacuum c. The time-dependent radius of the blackbody core Rc can be defined by the 



expression 
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R 

tot 



K'°\r,t)dr = 1, (35) 



Rc 



where is the wavelength-averaged absorption coefficient between 200 and 800 nm. Equa- 
tion (jHHjl implies that if radiation from a spherical surface at depth Rc is damped to some 
extent (the optical depth being 1), then radiation from further interior is completely opaque 
to an outside observer. The spherical surface Rc serves as the surface of a blackbody in place 
of radiations from the interior. The determination of Rc starts from the outermost. If 
is sufficiently large, there will be a finite-size blackbody core < Rc < R such that Eq. (j35p 
is satisfied; if is small, the left hand side of Eq. will be less than 1 even if Rc=0, 
implying that the bubble is optically thin. The calculated photon absorption coefficients 
indicate that the bubble is always optically thin, Rc = 0. However, in order to see how a 
finite-size blackbody model behaves, we intentionally amplify k*"* by a free parameter Ec so 
that 

/•R 

Ec-R'°\r,t)dr = 1 (36) 



Rc 

will give a nonzero Rc during the collapse stage. It is evident that larger Ec makes Rc 
closer to R. With Rc{t) at hand, we can calculate the light emission by Eq. (j33|) together 
with original whose small quantity makes the second term in Eq. (jHHI) dominant. For 
convenience of discussion we denote -Rc = as the optically thin model and Rc > as the 
finite-size blackbody model. It is remarked that the finite-size blackbody model is physical 
ii Ec = 1, and is ad hoc if Ec > 1. For fitting purpose, Ec is different from case to case but 
remains fixed during an acoustic cycle. 

It is meaningful to look at the light pulses and spectra. The integration of the spectral 
radiance over a suitable wavelength intervals (Auv = 200 nm ~ = 800 nm) gives the total 
power emitted into the measurable part of the spectrum, and integration over one acoustic 
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period Tg gives the light spectrum that can be compared with the experimental results 

p^\t) = f ^ pi\t)d\ = 1 r prmt (37) 

-'Auv Jo 
IV. NUMERICAL RESULTS 

It is well known that the ambient bubble radius Rq depends on experimentally controllable 
parameters such as the driving pressure amplitude Pa, the water temperature Too, and the gas 
concentration dissolved in the water Cqo- A problem with past SL spectrum measurements 
was they seldom gave the key parameters Pa and Rq at the same time. This left freedom 
for theoretical studies to fit experimental data using different Pa and Rq. The present study 
tries to use the same parameters as those in previous literatures. For comparison with 
other calculations, the parameters used here are Pa = 1.2 bar. Too = 298 K (Storey [3|) 
and Pa = 1.35 atm. Too = 300 K (Xu for identical equilibrium radius Rq = 4.5 /im 
and driving frequency f = 26.5 kHz, and Pa = 1.4 bar for Rq = 6.0 /xm, / = 20.6 kHz, 
Too = 293.15 K ( Moss [12^). The dissolved gas concentration is Coo/cq = Pno^/ Pq = 0.395% 
(3 Torr partial pressure) for all above cases. For comparison with the experiment Rq = 
4.5 fim (He), 5.5 fim (Xe), / = 42 kHz, Too = 296.15 K, P^^^ = 150 Torr (He) and 3 Torr 
(Xe), while Pa is adjustable. Other parameters are Poo = 101325 Pa, pioo = 996.6 kgm~^, 
ki = 0.609 Wm-^K-i, Cpi = 4179 Jkg"^ R-^, and A = 2 x lO'^ mVs. The van der 
Waals excluded volumes are given in Table I. Initial number densities of ions and electrons 
a.e e.t„„ated using the Saha equation Q. The ...tial bubble content eontain. 2% n,o.a. 
fraction water vapor. This number seems arbitrary, but our results are based on the second 
acoustic cycle when initial disturbances are presumed to be decayed. 

TABLE I: van der Waals excluded volumes. 



Species 


He & ions 


Ar &: ions 


Xe & ions 


H2O 


OH 


H 


H^ 





0+ 


H2 


O2 e- 


bi (cm'^/mol) 


23.7 


32.2 


51.0 


30.5 


15.25 


4.98 


36. 


8 2.77 


27.7 


26.6 


31.8 0.0 
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A. Effects of chemical reactions 



Figure 1 shows one forcing period of the radius of an argon bubble corresponding to 
case I (pure noble gas in the bubble), case II (with water vapor but no chemical reactions), 
and case III (with water vapor and chemical reactions) as labeled in Ref. j^]. One can 
see that the difference between II and III is almost indiscernible, but that between I and 
the latter two is large. The existence of water vapor increases the maximum radius and 
delays the collapse. In spite of little difference between II and III in the R-t curve, large 
difference occurs for thermodynamic variables at collapse. Table II shows comparison of 
some qua ntities. The maximum radius of the present calculation is smaller than that of 
Story [26[. This is mainly because the acoustic forcing terms in Eq. (I23|) is separate from 



Hh rather than absorbed in Poo as treated in Refs. [26|, |35| that magnified Pa by a factor 
of n/{n — 1) [see Eq. ()25|) ]. Since there are much differences between the present model 
and Storey's, quantitative discrepancies are expectable for the extreme values. Both models 
predicted temperature reduction from case I to III. However, the present result show that 
temperature is slightly reduced from I to II, but heavily from II to III. The slight reduction 
is due to the compensating effects of increased compression ratio (-Rmax/-Rmin) and reduced 
ratio of specific heats Q|. The larger reduction from 17 000 K (case II) to 8900 K (case III) 
indicates the significant effect of chemical reactions in reducing the temperature. 



TABLE II: Comparison of extreme values for an argon bubble Rq = 4.5 /xm, Pa = 1.2 bar. The 
amount of water vapor in II is in mole fraction and evaluated at the moment of -Rmin- 





Rma,x/ Pmin 




^max (K) 




total vapor (%) 


case 


present 


Ref. [26] 


present 


m 


present [26] 


I 


25.4/0.88 


28.0/0.80 


17900 


20900 




II 


28.7/0.76 


31.3/0.70 


17000 


9700 


7.7 W 


III 


28.9/0.70 


31.7/0.65 


8900 


7000 






2^ did not specify the exact moment for this value. 



The effects of chemical reactions on thermodynamic variables are best reflected in the 
distributions of temperature and chemical yields inside the bubble. Figure 2 shows snapshots 
of the spatial profiles of temperature around the moment of minimum bubble radius. It can 
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be seen that temperatures are considerably reduced in the reacting case III. Figure 3 shows 
the numbers of molecules of different species and temperature at the bubble center as a 
function of time. Note that in the first acoustic cycle the water vapor begins to dissociate 
appreciably at t = —2 ~ —1 ns, while in the beginning of the second cycle there are already 
some amounts of H2 and O2 gases accumulated. The chemical reactions occur in a time 
scale of several nanoseconds, producing considerable amounts of H, O, OH radicals and H2 
and O2 gases. It is remarked that the prediction of chemical products is very difficult as 
the reaction mechanisms and phase change processes are largely unknown under extreme 
conditions in a SL bubble. 

Next we compare a -Rq = 6.0 fim Ar bubble driven at Pa = 1.4 bar, which was labeled as 



Al in Ref. 



12| . The present temperatures at the bubble center are respectively 109 600 K, 



34 700 K, and 16 200 K for cases I— III in Fig. 1, suggesting that both the reduced ratio 
of specific heats due to the presence of water vapor and the chemical reactions contribute 
significantly to the reduction of temperatures. The amount of trapped water vapor at the 
moment of -Rmin for case II occupies 23% molar fraction, smaller than 33 % [2^. Figure 4 
shows several snapshots of the spatial profiles of thermodynamic variables. A main feature 
is that only compression waves occur. As seen from the velocity profile, a compression 
wave moves outward at and ts, refiects from the bubble wall and moves inward at t^. 

waves were reported. A 

possible reason is that the formation of a shock is sensitive to differences in equations of 
state, accommodation coefficients, chemical reactions, and treatments of the liquid motion, 
and so on. Another feature in Fig. 4 is that the temperature in the inner zone is reduced 
more severely than in the outer zone at ti,t2,ti, and tg when the compression wave is not 
strong. This is because considerable water vapor is trapped in the inner zone of an Ar bubble 



result of thermal diffusion 
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|35| (which states that a heavier species tends to diffuse 
toward the cooler region) and dissociates there, thus peak temperature is not at the center, 
but at some place close to the bubble interface. For a lighter-than-water-vapor He bubble, 
less water vapor is congregated in the central zone, and temperature peak will be located 
at the bubble center as [will be shown in Fig. 6 (b)]. Figure 5 shows one snapshot of the 
number density distributions and the degrees of ionization. It can be seen that the amounts 
of products due to chemical reactions are considerable, while the degrees of ionization are 
quite small (the maximum being 2.8% for O^), contrary to significant ionizations when water 
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vapor was not taken into account jl5|. The degrees of ionization of H and O far exceed that 
of Ar although Ar atom is more ample in quantity. 

The numerical results in this subsection demonstrates that the trapped water vapor and 
the ensuing endothermic chemical reactions significantly reduce the temperature, resulting 
in very low degrees of ionization. The chemical reactions can produce considerable amounts 
of atomic and molecular species, some of which, such as H and O atoms, are easier to ionize 
than a noble gas atom such as He or Ar. As the evaporation is a robust process in the 
bubble oscillation, chemical radicals will have significant influence on SBSL mechanism as 
will be shown in Sec. IV C. 



B. Effects of noble gas types 

Previous numerical studies pointed that shock formation depends sensitively on, among 

n nn n n 

other factors the amount of water vapor and its distribution |2a|. Xu et al. |31| 

showed that shock waves develop in a bubble filled with 70% Xe and 30% (mole fractions) 
water vapor, but no shocks occur for similarly filled Ar or He bubbles. With evaporation- 
condensation process and chemical reactions taken into account, we are able to investigate 
effects of noble gases on the thermodynamic processes more realistically. We calculated Xe, 
Ar and He bubbles using the same i?o, -Pa, /, and Too as those in Ref. |3ll |. 

Figures 6(a) and 6(b) show snapshots of the spatial profiles of thermodynamic variables 
around the moment of minimum bubble radius for Xe and He bubbles. Snapshots of the Ar 
bubble are similar to those in Fig. 4, thus are not shown here. In the Xe bubble [Fig. 6(a)], 
it is seen that an inward-going compression wave at t2 evolves into a strong outward-going 
shock at t^. The first focusing of the shock happened between t2 and ^3 leads to extreme high 
temperatures (> 10® K), but the duration is very short (< 1 ps) and the region is confined 
to the center (r < 0.005 /im). However, Fig. 6(b) shows that only weak compression waves 
occur in the He bubble. These results are in qualitative agreement with those of Xu et 
al. [3^. Note that the temperature peaks are often at the center in He bubble except when 
a wavy disturbance reflects from the bubble wall at t^. This feature mainly results from 
the thermal diffusion between the light He gas and the heavy water vapor as mentioned in 
previous subsection. 

Figure 7(a) and 7(b) show one snapshot of the number density and the compositional 
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distributions at the moment of minimum bubble radius for Xe and He bubbles. In both 
bubbles there are significant numbers of chemical products, especially H, O, OH, and H2. 
However, the right figure in Fig. 7(a) indicates that atomic species in Xe bubble are sig- 
nificantly ionized only in the central zone (r < 0.07i?), while the right figure in Fig. 7(b) 
shows that the degrees of ionization in He bubble are small, especially that for He gas. The 
situation for Ar bubble is found to be in between [Fig. 5(b)]. This suggests that chemical 
products of the water vapor may play an increasingly important role in SBSL from Xe, Ar, 
to He bubbles. 



C. Calculated light spectra and pulses 

We shall compute the emitted lights by using the optically thin model and the ad hoc 
inite-size blackbody model. We fit our calculations to Fig. 2 of a recalibrated experiment 
Ifij ] under the same parameters as given in the beginning of this section. The former model 
has only Pa while the latter has Pa and Ec in Eq. as fitting parameters. 

Figure 8 shows comparison of the spectral radiances. It is seen that the optically thin 
model does not match well with the experimental spectrum of Xe bubble, but the finite-size 
blackbody model matches well. The fitting Pa = 1.55 atm of the optically thin model seems 
to be out of the stable SBSL range (1.2 — 1.5 atm), while the fitting Pa = 1.28 atm of the 
ad hoc finite-size blackbody is within the range, of course with the help of a large value 
Ec = 1.8 X 10^. (Fitting using Pa = 1.49 atm and a smaller value, E^ = 60 can also give 
similar spectrum but the resulting FWHM is only 40 ps, much shorter than experimental 
200 ps.) However, it can be seen that either models are unable to fit the spectrum of He 
bubble. The finite-size blackbody is much better than the optically thin model as the latter 
deviates severely from the experiment. The maximum temperatures at Pa = 1.28 atm for 
Xe bubble and at Pa = 1.45 atm for He bubble are 8600 and 16 700 K, respectively, which 
are c omp arable to the blackbody fitting temperatures 8000 (Xe) and 20 400 K (He) used in 



Ref, 

Figure 9 compares the time variations of the normalized power for the "measurable," 
"UV" (300—400 nm), and "red" (590—650 nm) wavelength intervals [4]. Both models show 
good wave-length independence of light pulse, a key ingredient of SBSL thought by several 
researchers 0,1^. One curves "A," the optically thin model gives FWHM of 26 ps for Xe 
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bubble, and 13 ps for He bubble, which are much shorter than the experimental flash widths 

M 

of 200 ps (Xe) and 100 ps (He) 116]. On curves "B," the finite-size blackbody model predicts 
FWHM of 185 ps for Xe bubble, and 18 ps for He bubble, better than the optically thin 
model. Figure 10 shows the variations of the blackbody core and bubble radius with time. 
One can see that the blackbody core appears abruptly, attains maximum around the moment 
of minimum bubble radius, and disappears suddenly. The short life of the blackbody core 
explains why this blackbody model also shows wave-length independence of the light pulse 
similar to the optically thin model: The quick rise and fall of Rc in accordance with the 
variation of the photon absorption coefficients cut down the long fall time for red light j^. 

Figure 11 shows the visible light powers contributed from the total and partial species 
as computed by the optically thin model. In Fig. 11(a) we see that the power from xenon 
is dominant, while that from the water vapor and its chemical products contributes a little. 
But in Fig. 11(b), the power from helium is small, while the water vapor and its chemical 
products contribute dominantly with H and O radicals being the primary ones. Although 
the absolute values of the light powers from the optically thin model are quite small, the 
relative contributions to the total power verify previous postulation that light emission from 
radicals generated from water vapor dissociation may dominate SBSL in the He bubble j^. 

The failure of the optically thin model to match with the experiment and the improvement 
by adding an ad hoc finite-size blackbody core to it raise contradiction. The major reason for 
the failure of the optically thin model is the significant reduction of temperature due to the 
existence of water vapor and endothermic chemical reactions. We remark that the present ad 
hoc blackbody is not the unique way for better fitting, but is shown better than other ad hoc 
ways such as multiplying the photon absorption coefficient by an arbitrary factor. We did 
not exclude the possibility that the optically thin model could fit well to the experimental 
data should higher temperature or increased opacity be obtained in whatever a natural 
way. As a numerical study, we have made tests of model sensitivity to parameters. The 
evaporation-condensation process and chemical reaction rates are important for the modeled 
;emperature. We used the accommodation coefficient au = 0.4 [26'] as well as the formula 



2J|, and found that the former led to higher temperature so that we could use Pa = 1.7 atm 



to obtain the same spectrum of He bubble as the one obtained by using Pa = 2.0 atm in 
Fig. 8. However, the results were still unable to fit to the data. As to the chemical reaction 
.ate. we eo.pa..ed two dffle.e. .eae«o. .ate sets QQ. TKe.e was sH,Kt difference . 
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quantity and what we presented here was the one corresponding to higher temperature js^ . 
We have aheady adopted the modified chemical equihbrium constant p2|) for a van der 
Waals gas. This modification was claimed to let the optically thin model give sufficient light 
emission {2^, but we found while it raised temperature to some extent, it could not result 
in enough light emission. Therefore, it will be better to pursue other ways to reconcile the 
conffict. One conjecture was if there existed a mechanism that would greatly increase the 
photon absorption coefficient of the highly compressed bubble content js], |l6| , such as the 
lowering of ionization potentials . it is also wished to have better theories to compute the 
photon absorption coefficients of a very dense gas mixture, and to take into account other 
light emission processes due to the existence of chemical products. Meanwhile further efforts 
are necessary to reduce modelling uncertainties such as the chemical reaction rates under 
high-pressure and high-temperature conditions, and the description of the surrounding liquid 
motion. The RP equation approach in this and other studies 12611 should be more critically 
compared with the full hydrodynamic equation approach 3, ll<^ . 



V. CONCLUSIONS 



A reffned hydrochemical model is presented to simulate the complex processes inside a 
sonoluminescing bubble. The numerical simulations of Xe, Ar, and He bubbles indicate that 
the trapped water vapor and its endothermic reactions reduce the temperature significantly. 
In the stable SBSL range, at most compression waves can appear in He or Ar bubbles, 
while shock waves can occur in Xe bubbles only for higher driving amplitudes. The lower 
temperature in the bubble rarely leads to appreciable ionization except for Xe bubble at 
the center during the shock wave focusing. The chemical radicals generated from water 
vapor dissociations become increasingly important in the light emission from Xe, Ar, to He 
bubbles. Particularly, H and O radicals are shown to be the primary light-emitting matters 
in He bubbles. 

The key finding of this study is that the optically thin thermal emission model was unable 
to match with experimental data mainly due to the reduced temperatures in the bubble. 
The introduction of a finite-size blackbody core made the calculated light spectra and pulse 
widths match better with experimental ones. The present expertise to define an optically 
thick region is ad hoc and only serves to illustrate one possible improving direction. Im- 
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provements by physical considerations such as the lowering of ionization potentials and the 
refinement of the photon absorption processes under the extreme conditions of sonolumines- 
cence, and by reduction of modeling uncertainties, are worthy of further investigation. 
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FIG. 1: Radius of an argon bubble vs time over one acoustic period for the same parameters as used 
in Ref. 2g]: Rq = 4.5 /xm. Pa = 1.2 bar, / = 26.5 kHz, T^o = 298 K. Case I: without phase change 



and chemistry; case II: with phase change but without chemistry; case III: with phase change and 
chemistry. 
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FIG. 2: Snapshots of the temperature distributions for Ar bubble at Pa = 1.2 bar, Rq = 4.5 //m. 
The upper three hues arc for nonrcacting case II, and the lower three lines arc for reacting case 
III. t = ps corresponds to the time of minimum radius (imin = 19.015945 /xs). 
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FIG. 3: The numbers of molecules of species and the temperature of the center around the moment 
of minimum bubble radius as a function of time for Ar bubble at Pa = 1.2 bar, Rq = 4.5 /xm. (a) 
First acoustic cycle, (b) second acoustic cycle, (c) temperature at the bubble center. 
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FIG. 4: Snapshots of the spatial profiles of density, velocity, pressure, and temperature for Pa = 
1.4 bar, Rq = 6 /xm. Time sequences are ti = —100 ps, t2 = —80 ps, ts = —60 ps, t4 = —40 ps, 
= —20 ps, = ps, h = 40 ps, tg = 100 ps, where = denotes the time of minimum radius 
{tmin = 28.970309 /is). 



26 



E 

o 
E 



"a 
o 

E 
3 



(a) 



1 1 1 1 1 1 1 


- 


























1 , 1 





0.2 



0.4 
r(|xm) 



0.6 



0.8 




FIG. 5: The spatial profiles of number densities for molecular species (a) and degrees of ionization 

for ionic species (b) at the time of minimum bubble radius t = (imin = 28.970309 /is) for 



Pa = 1.4 bar, i?o = 6 /xm. The degree of ionization is computed using the Saha equation. 
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FIG. 6: Snapshots of the spatial profiles of density, velocity, pressure and temperature for Pq = 
1.35 atm, Ro = 4.5 /xm. (a) Xe bubble. Time sequences are ti = —80 ps, t2 = —40 ps, = —20 ps, 
t4 = ps, t5 = 40 ps, te = 60 ps, where t4 = denotes the time of minimum radius (tmin = 
22.006191 /Its), (b) He bubble, h = -150 ps, t2 = -100 ps, ts = -50 ps, i4 = ps, t5 = 100 ps, 
te = 250 ps, where tmin = 22.057459 fis. 
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FIG. 7: The spatial profiles of the number densities of molecular species (left) and the compositions 
of ions (right) for Pa = 1.35 atm, Rq = 4.5 /xm. (a) Xe bubble with ions computed using the 
nonequilibrium ionization, (b) He bubble with the degree of ionization computed using the Saha 
equation. All species are shown for the time of minimum radius. 
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FIG. 8: Spectral radiance of the SL light from bubbles of Xe and He in water. The same parameters 
as in experiment IQ] are Xe (the ambient radius Rq = 5.5 fim, dissolved partial pressure 3 Torr) 
and He (i?o = 4.5 fiui, 150 Torr), water temperature 23 °C, driving frequency 42 kHz. The squares 
and triangles are experimental spectra of Xe and He bubbles, respectively. The solid lines are 
calculated spectra of the ad hoc finite-size blackbody model with fitting parameters Pa = 1.28 atm, 
Ec = 1.80 X 10^ (Xe), and Pa = 1.45 atm, Ec = 2.50 x 10^ (He), and the dashed lines are those of 
the optically thin model with fitting parameters Pa = 1.55 atm (Xe) and 2.0 atm (He). 
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FIG. 9: Normalized radiation power vs time for the optically thin model ("A") and the ad hoc 
finite-size blackbody model ("B") in Xc bubble (a) and He bubble (b). The solid line denotes the 
total measurable power, the dashed line in the UV range (300 nm < A < 400 nm), and the dotted 
line in the red range (590 nm < A < 650 nm) . Time is relative to the moment of minimum bubble 
radius. The parameters are the same as in Fig. 8. 
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FIG. 10: The blackbody core Rc (solid line) and the bubble radius R (dashed line) vs time for Xe 
and He bubbles. The parameters are the same as those for the ad hoc finite-size blackbody model 
in Fig. 8. 
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FIG. 11: The optical powers vs time for Xe (a) and He (b) bubbles computed from the optically 
thin model. The uppermost curve is the total power, and others are powers contributed from the 
species as marked. The driving pressure amplitudes are 1.28 atm for Xe bubble, and 1.45 atm for 
He bubble. Other parameters are the same as those in Fig. 8. Time is relative to the moment of 
minimum bubble radius. 
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